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Mid-Infrared Source Multiplicity within Hot Molecular Cores 
traced by Methanol Masers 



o 
o 

(N 

Oh: 

< 

(N 



> 

(N 

o 
o 

o 



S. N. Longmore^'^, M. G. Burton\ V. Minier^ and A.J.Walsh^ 

^School of Physics, University of New South Wales, Sydney, 2052, NSW, Australia 
^Australia Telescope National Facility, CSIRO, Epping, 1710, Sydney, Australia 
^Sen-ice dAstrophysique, DAPNIA/DSM/CEA Saclay 91191 Gif-sur-Yvette, France 



Accepted 2006 March 22. Received 2006 February 6 



ABSTRACT 

We present high resolution, mid-infrared images toward three hot molecular cores signposted 
by methanol maser emission; G173.49+2.42 (S231, S233IR), G188.95+0.89 (S252, AFGL- 
5180) and G192.60-0.05 (S255IR). Each of the cores was tai-geted with Michelle on Gemini 
North using 5 filters from 7.9 to 18.5 /im. We find each contains both large regions of ex- 
tended emission and multiple, luminous point sources which, from their extremely red colours 
(^18.5/^7.9 ^ 3), appear to be embedded young stellar objects. The closest angular separa- 
tions of the point sources in the three regions are 0.79, 1 .00 and 3.33" corresponding to linear 
separations of 1,700, 1,800 and 6,000 AU respectively. The methanol maser emission is found 
closest to the brightest MIR point source (within the assumed 1" pointing accuracy). Mass 
and luminosity estimates for the sources range from 3-22 Mq and 50-40,000 L0. Assuming 
the MIR sources are embedded objects and the observed gas mass provides the bulk of the 
reservoir from which the stars formed, it is difficult to generate the observed distributions for 
the most massive cluster members from the gas in the cores using a standard form of the IMF. 

Key words: masers - stars:formation - techniques: high angular resolution - stars:early type 
- stars:mass function - infrared:stars. 
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1 INTRODUCTION 

Massive stars play a fundamental role in driving the energy flow 
and material cycles that influence the physical and chemical evolu- 
tion of galaxies. Despite receiving much attention, their formation 
process remains enigmatic. Observationally, the large distances to 
the nearest examples and the clustered mode of formation make it 
difficult to isolate individual protostars for study. It is still not cer- 
tain, for instance, whether massive stars form via accretion (similar 
to low mass stars) or through mergers of intermediate mass stars. 

Advances in instrumentation, have enabled (sub) arcsecond 
resolution imaging at wavelengths less affected by the large col- 
umn densities of material that obscure the regions at shorter wave- 
lengths. Recent observations exploiting these capabilities have un- 
covered the environment surrounding individual massive protostel- 
lar system s . Fr om analysis of ~2.3 /im CO bandhead emission, 
lBik& Thi (2004) have inferred Keplerian disks very closely sur- 
rounding (within a few AU) four massive young stellar objects, 
while interferometric, mm-continuum observations, find the mass- 
function of protostellar dust clumps lies close to a Salpeter value 
down to clump radii of 2000AU (Beuther & Schilke 2004). These 
high resolution observations point toward an accretion formation 
scenario for massive stars. 

Further discrimination between the two competing models is 
possible by examining the properties, in particular the young stel- 



lar populations, of hot molecular cores. The mid-infrared (MIR) 
window (7-25 /^m) offers a powerful view of these regions. The 
large column densities of material process the stellar light to in- 
frared wavelengths, and diffraction limited observations are readily 
obtained. 



Recent observations indicate that class II met hanol masers 
exclusively trace regions of massive star formation JMinier et alJ 
2003 ) and are generally eith er not associa t ed or o ffset from UCHii 
regions iWalsh e t al. 1998). iMinier et alj pOO?) (hereafter M05) 
have carried out multi-wavelength (mm to MIR) observations to- 
ward five star forming complexes traced by methanol maser emis- 
sion to determine their large scale properties. They found that 
maser sites with weak (<10mJy) radio continuum flux are associ- 
ated with massive (>50Mq), luminous (>1O*L0) and deeply em- 
bedded (Au >40 mag) cores characterising protoclusters of young 
massive (proto)stars in an earlier evolutionary stage than UCHii re- 
gions. The spatial resolution of the observations (>8") was, how- 
ever, too low to resolve the sources inside the clumps. Details of 
the regions from observations in the literature are described in M05. 
We have since observed three of the M05 regions at high spatial res- 
olution to uncover the embedded sources inside the cores at MIR 
wavelengths. 
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2 OBSERVATIONS AND DATA REDUCTION 

The data were obtained with Michelle^ on the 8-m, Gemini North 
telescope in queue mode, on the 18* \ 22" and 30* of March 
2003. Each pointing centre was imaged with four N band silicate 
filters (centred on 7.9, 8.8, 11.6 and 12.5 /im) and the Qa filter 
(centred on 18.5 ^m) with 300 seconds on-source integration time. 
G173.49 and G188.95 were observed twice on separate nights and 
G192.60 observed once. 

The N and Q band observations were scheduled separately due 
to the more stringent weather requirements at Q band. The standard 
Chop-Nod technique was used with a chop throw of 15" and chop 
direction selected from MSX images of the region, to minimise off- 
field contamination. The spatial resolution calculated from standard 
star observations was ~ 0.36" at 10 /^m and ~ 0.57" at 18.5 /im. 
The 32"x24" field of view fully covered the dust emission observed 
by M05 in each region. 

Particular care was taken to determine the telescope pointing 
position but absolute positions were determined by comparing the 
MIR data to sensitive, high resolution, cm continuum, VLA images 
of the 3 regions (Minier et al. in prep). Similar spatial distribution 
and morphology of the multiple components allowed good registra- 
tion between the images. The astrometric uncertainty in the VLA 
images is ^^1". 

Flux calibration was performed using standard stars within 
0.3 airmass of the science targets. There was no overall trend in 
the calibration factor as a result of changes in airmass throughout 
the observations. The standard deviation in the flux of standards 
throughout the observations was found to be 7.4, 3.1, 4.4, 2.4 and 
9% for the four N-band and 18.5 ^m filters respectively. The statis- 
tical error in the photometry was dominated by fluctuations in the 
sky background. Upper flux limits were calculated from the stan- 
dard deviation of the sky background for each filter and a 3(t up- 
per detection limit is used in Table 1. Similarly, a 3o- error value 
is quoted for the fluxes in Table 1 (typical values for the N and Q 
band filters were 0.005 and 0.03 Jy respectively). The flux densities 
for the standard stars were taken from values derived on the Gemini 
South instrument, T-ReCS^ which shares a common filter set with 
Michelle. 

Regions confused with many bright sources were deconvolved 
using the Lucy-Richardson algorithm with 20 iterations. This was 
necessary to resolve source structure and extract individual source 
fluxes. The instrumental PSF was obtained for each filter using a 
bright, non-saturated standard star. The results were reliable and re- 
peatable near the brighter sources when using different stars for the 
PSF and observations of the objects taken over different nights. As 
a further check, the standard stars were used to deconvolve other 
standards and reproduced point sources down to 1% of the peak 
value after 20 iterations, so only sources greater than 3% of the 
peak value were included in the final images. The resulting decon- 
volutions are shown in Fig 1. 

3 DERIVING PHYSICAL PARAMETERS 

A colour temperature was derived for each object, assuming the 
flux was greybody emission from dust at a single temperature 
and emissivity given by 1 — e~^^. tx was calculated using Ai, 
(from M05) and the IPraind (l989|) normalised silicate profile 



^ The Michelle imager employs a Raytheon 320x240 pixel Si:As IBC array 
with a pixel scale of 0.099" . 
^ See http ://w w w. gemini.edu | 



(T"i2Mm=Ti8Mm.=0.022^u). The 12.5 to 18.5 /im flux ratios were 
used to derive the temperature as they are expected to contain the 
least PAH emission and are less affected by the broad-band silicate 
absorption features. In cases where either the 12.5 and 18.5 ^m 
source flux was below the detection limit, only an upper and lower 
temperature limit respectively was derived. Luminosities were de- 
rived in the blackbody limit at the calculated source colour temper- 
ature. The values agreed to <1% with those calculated by integrat- 
ing the product of the Planck function and the emissivity function 
from 1-600/xm so provide a reasonabl e lower limit to the bolomet- 
ric luminosity JDe Buizer et alj|2005l) . Luminosities derived from 
colour temperatures agree well with those from radiative trans- 
fer model fits to the SED across the near and MIR (De Buizer, 
priv. comm.). A mass estimate for the sources was calculated by 
considering the range of plausible a values in a mass luminosity- 
relationship, LocM", for a deeply-embedded, zero age main se- 
quence sta r. For stars on the main sequence <2OM0, q=3.45 
(eg lAUerJ (1975)), while for stars > 20Mq, a ~3 JZinneckeil 
I2OO3I) . Additional luminosi ties due to accretio n (accounting for up 
to half the total luminosity iOsorio et alll999l) ) may increase a up 
to ~4. Over this range of a (3<a<4), a factor of ten error in lu- 
minosity corresponds to only a factor of ~2 in the calculated mass. 
An intermediate value, a=3.45, was chosen to derive the mass es- 
timate. Flux errors of 20% and a 20% distance uncertainty corre- 
sponds to a colour temperature error of 10% and luminosity error 
of 35%. However, the relative luminosity uncertainty for sources 
within a core drops to 20% as they lie at the same distance. Using 
the 3a upper flux limits from the 12.5 and 18.5/im images the lower 
mass detection limit is ~ IM0 for an object at a typical observed 
source temperature of 130K at 1.8 kpc. The source positions, fluxes 
and derived properties are shown in Table 1. 



4 RESULTS AND DISCUSSION 

Figure 1 shows the 18.5 and 7.9 /im images of the three regions. 
The insets show deconvolved 7.9 jim images toward bright, con- 
fused regions. All the sources toward each region are assumed to 
belong to the same protocluster due to their sharply rising MIR 
flux (Fis.s/Fr.g > 3) and large optical depths (yl„ > 40 mag). 
Although F18.5 ~ Fj.g for G192.60:2 1, the spectral energy dis- 
tribution peaks in the MIR so it must still be heavily embedded. 
Due to the low signal to noise of the detection G 173. 49: 1 5 was not 
included in the following analysis. 

Table 1 lists the positions, morphology, measured fluxes 
and derived properties of the detected MIR sources and Table 2 
shows the large scale clump properties derived in M05. In all three 
clumps, the individual source radii derived from the black body 
assumption compare well to the observed radii. The variation in 
individual source radii and temperature compared to the M05 
values in Table 2 is not unexpected due to the M05 assumption 
of a single powering so urce. With a factor of ~2 error in the 
M05 luminosity (see eg IPurcell et alJ i200Q for error analysis 
of two-component, greybody fitting) the total measured MIR 
luminosity (Lmir) for G188.95 and G192.60 accounts for the M05 
total core luminosity (Ltot)- The large discrepancy in G173.49, 
however, suggests there are undetected, luminous, cold source(s). 
Attributing the missing flux to a single source would correspond to 
a star of mass '^20Mq . 
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0.9±0.06 


<51 


- 


- 


- 



Table 1. Positions, measured fluxes and derived properties of the detected mid-infrai'ed objects in the three cores. T^ol is th^ colour temperature of the 
source derived from the 12.5 and 18.5 /xm flux ratios. R and Lmi,. are the radius and luminosity of a greybody at temperature T^oi- The table uses notation 
Oxxx.xx:a b for the sources, where Oxxx.xx is the formation region from M05, 'a' is the clump number from M05 and 'b' is the mid-infrared source number 
The superscript p, pe and e in the first column, denote point sources, point sources with low level surrounding extended emission and extended sources 
respectively. The superscript m denotes the source closest in projection to the methanol maser. The radius in the final column is from the black body fit and 
is consistent with observations. The quoted errors and upper flux limits are 3(t values. Distances of 1.8, 2.2 and 2.6 kpc are used for 0173.49, G188.95 and 
0192.60 respectively. 



4.1 Nature of the MIR sources 

MIR emission traces dust heated by nearby sources. Its morphology 
provides clues as to the nature of the source, in particular, whether 
the dust is heated externally or by an embedded source. Observa- 
tionally, the sources fall into three morphological categories; (i) un- 
resolved point sources, (ii) unresolved point sources with weak sur- 
rounding extended emission and (iii) extended emission, which we 
have designated in Table 1 as P, PE and E, respectively. The source 
morphology is consistent across the filter range (after deconvolu- 
tion at longer wavelengths) but it is unclear whether G 173. 49: 1 3 is 
actually a colder component of the extended emission surrounding 
G173.49:l 2 instead of a separate source. 

With a spatial resolution ~400AU we have made the assump- 
tion that the point-like sources (P and PE) are stellar in origin and 
the extended sources (E) are non-stellar (eg externally heated) with 
the exception that the most luminous source in the region is ex- 
pected to be internally heated regardless of morphology. The inter- 
nal heat sources are expected to be either very young stars or proto- 
stars. AU non-stellar sources have been excluded from the analysis 
in §§4.2, 4.3 and 4.4. 



4.2 Point source spatial distribution 

As the large scale clump dust and gas morphology appears simple 
and centrally peaked (see M05), we make the reasonable assump- 
tion that the protocluster centres coincide with the central peak of 
dust emission. The spatial distribution of the point sources within 
the protocluster is similar between the clumps with close point 
sources toward the cluster centre. The methanol masers are found 
closest to the brightest MIR point source (within the assumed 1" 
pointing error from image registration). These sources have tem- 
peratures sufficient to evaporate methanol ice from the dust grains 
into the gas phase (>90K) as well as sufficient luminosity of IR 
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Table 2. Properties of the three cores derived from M05 using a greybody 
approximation to fit the spectral energy distribution. Ltot is the total clump 
luminosity, L^oj, T^ot 2nd R^ot give the luminosity, temperature and ra- 
dius of the hot component derived from the greybody approximation. h,„ir 
is the total luminosity of the sources in each core from Table 1. Mgas is 
the average gas mass derived from mm-continuum and molecular spectral 
line observations of the clumps. The NH2 column densities are the values 
derived from mm-continuum observations in M05. 



photons to pump the masing transition - conditions m odels suggest 
are required for such emission tVoronkov et all2005l) . 

It is known that mo re massive stars favour cluster centres 
(e.g. 'de G riis et alj i2002h ). but it is unclear whether they form 
there or migrate in from outside. We have used the simple-harmonic 
model of ballistic motion developed bv lWalsh et alj ^200^ to con- 
sider the motion of sources within the cores. Using the measured 
column density and radius from M05 (listed in Table 2), the time 
required for migration from the edge to the centre is ~ 10^ years. 
This is comparable to the predicted HMC lifetime of 10^ years 
(Kurtz et al. 2000) so we can not rule out the possibility of migra- 
tion within the clumps. Any sources having migrated to the centre 
in this way would have acquired a velocity of '--^ 2 kms~^ with re- 
spect to the clump. 
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4.3 MIR source multiplicity 

Massive star s in clusters are observed to have a h igh companion 
star fr action iZirmeckeill2003h . In the M16 cluster. IPuchene et alj 
i200l ll) observed massive stars (earlier than B3) with visual com- 
panions separated by 1000-3000AU. If multiple systems are bound 
from birth, it is likely some of the sources we have observed will 
belong to multiple systems, even though the companions may lie 
below the detection limit. However, all three regions show two or 
more point sources at close angular separation (see insets of Figure 
1) corresponding to linear separations of 1700 to 6000AU. We can- 
not determine whether these stars are physically bound or simply 
close due to projection effects but we can calculate the instrumen- 
tal sensitivity required to confirm or deny the association. Assum- 
ing they are physically bound in a Keplerian orbit, the maximum 
proper motions (projection angle = 0°) of ~ 0.1 mas/year are too 
small to be detected on short temporal baselines. The maximum ve- 
locity difference (projection angle = 90°) ^ 2 kms"'^ is achievable 
by high spectral resolution observations of any line features. 



• Multiple, MIR sources which can be separated into three mor- 
phological types: unresolved point sources (P); unresolved point 
source with weak surrounding extended emission (PE) and ex- 
tended sources (E). 

• The point sources lie at close angular separations. Future high 
spatial and spectral resolution observations may be able to deter- 
mine whether or not they are physically bound. 

• The methanol masers are found closest to the brightest MIR 
point source (within the assumed 1" pointing accuracy). 

• Cooler, extended sources dominate the luminosity. 

• The time scale for a source at the core edge to migrate to the 
centre is comparable to the Hot Molecular Core lifetime, so it is not 
possible to rule out large protostellar motions within the core. 

• From the derived gas mass of the core and mass estimates for 
the sources, Monte Carlo simulations show that it is difficult to gen- 
erate the observed distributions for the most massive cluster mem- 
bers from the gas in the core using a standard form of the IMF. This 
conclusion would not hold, however, if most of the original gas has 
already formed stars, or has been dispersed such that the original 
core mass is much greater than now observed. 



4.4 Protocluster mass distribution 

The mass distribution of stars is generally well described as a power 
law through the initial mass function (IMF). Given the mass of gas 
available to form stars, we may estimate the likelihood that a clus- 
ter will end up with the most massive stars that are observed in 
it. The fraction of gas that forms stars is given by the star forma- 
tion efficiency (SEE) and is observationally found to be less than 
50% in any cloud an d to be ^ 33% for nearby embedded clusters 
iLada & Lada l2003'). For a cluster whose total stellar mass is 120 
(50, 3 20) M(7) (e quivalent t o the g as mass determined for the three 
cores). IWeidner & K roupa 1 2004 ) estimate that the mean maximum 
mass that a star may have in it is 10 (5, 20) M0 . This is comparable 
to the largest observed mass in two out of the three cases. However, 
we also observed several other stars in each cluster so can estimate 
the probability of generating stars of equal or greater mass than 
the remaining mass distribution. We did this by running Monte- 
Carlo simulations to p opul ate 10^ cluste rs using Salpeter 1 1955i) . 
iMiller & Scalol J 19791) and lScald h98d) IMFs until the available 
gas mass was exhausted. We only considered clusters which con- 
tained a star of at least equal mass to the most massive observed. 
The simulations show that even using the Salpeter form of the IMF 
(most biased toward forming high-mass stars) and allowing 50% of 
the gas to form stars, it is difficult to generate the observed mass 
distributions (probabilities < 10~^, 10~^, 10~^ for the three cores 
respectively). By itself, this may not be significant for a single clus- 
ter. However, since the probability is low for all three sources stud- 
ied, it is unlikely that the mass distribution of the most massive stars 
can be produced by sampling a standard form of the IMF from the 
reservoir of gas available for star formation. This conclusion would 
not hold if there was a substantial stellar mass already in the cluster 
that remains unseen, or if much of the original gas mass had already 
been dispersed from the core due to star formation. The former re- 
quires a SEE close to unity and given the relatively quiescent state 
of the cores, the latter seems unlikely. 



5 CONCLUSIONS 

A larger sample of young, massive protoclusters is required to 
draw general conclusions. However, in all three hot molecular cores 
traced by methanol maser emission we have found: 
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Figure 1. Images of the three regions. The top row displays the full field of view at 18.5 fira toward G173.49+2.42 Clump 1 (left), G188.95+0.89 Clump 1 
(centre) and G192.60-0.05 Clump 2 (right). The centre of the mm dust peak emission from M05 is shown as an ellipse, methanol maser sites are shown as plus 
symbols and radio continuum sources as boxes. A further ellipse in the G188.95+0.89 images shows the position of the IRAS point source within the field 
but registration with a single MIR source is not possible due to the much larger IRAS beam size. The next three rows show close up images at 18.5 /^m (left) 
and 7.9 /im (centre) of the bright point sources close to the dust peak and methanol maser emission (G173.49+2.42 (second row), G188.95+0.89 (third row) 
and G192.60-0.05 (bottom row)). The inset on the right shows the result when the boxed region in the 7.9 ^m image is deconvolved with the instrument PSF. 
Scale bars in the deconvolved images give the linear separation of the sources assuming they lie at the same distance. 



